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Abstract. The well-studied VHE (E > 100 GeV) 
blazar Mrk 501 was observed between March and 
May 2008 as part of an extensive multiwavelength 
observation campaign including radio, optical, X-ray 
and VHE gamma-ray instruments. Mrk 501 was in 
a low state of activity during the campaign, with 
a low VHE flux of about 20% the Crab Nebula 
flux. Nevertheless, significant flux variations could be 
observed in X-rays as well as 7-rays. Overall Mrk 
501 showed increased variability when going from 
radio to 7-ray energies. 

The broadband spectral energy distribution during 
the two different emission states of the campaign 
was well described by a homogeneous one-zone 
synchrotron self-Compton model. The high emission 
state was satisfactorily modeled by increasing the 
amount of high energy electrons with respect to the 
low emission state. This parameterization is con- 
sistent with the energy-dependent variability trend 
observed during the campaign. 
Keywords: Blazar, Mrk 501, SSC model 

I. Introduction 

Mrk 501 is a well-studied nearby (redshift z = 0.034) 
blazar which was first detected at TeV energies by 
the Whipple collaboration in 1996 JT]. In subsequent 
years Mrk 501 was regularly observed and detected 
in VHE 7-rays by many other Cherenkov telescope 
experiments. In particular during the whole year 1997 

°For a full author list, please see R. Ong et al (these proceedings) 
or http://veritas.sao. arizona.edu/conferences/authors ?icrc2009 ' (VER- 
ITAS) and http://wwwmagic.mpp.mpg.de/collaboration/members/ 
(MAGIC) 



when it showed an exceptionally strong outburst with 
peak flux levels up to 10 times the Crab Nebula flux 
and flux-doubling time scales down to 0.5 days J2j- 
Mrk 501 also showed strong flaring activity at X-ray 
energies during that year. The X-ray spectrum obtained 
was very hard and the synchrotron peak was found to be 
at ~ 100 keV, about 2 orders of magnitude higher than 
in previous observations 0. In the following years, Mrk 
501 showed only low 7-ray emission (of the order of 20- 
30% the Crab Nebula flux), apart from a few single flares 
of higher intensity. In 2005, the MAGIC telescope was 
able to observe Mrk 501 during another high-emission 
state which, although at a lower flux level compared to 
1997, showed flux variations of an order of magnitude 
and unprecedented flux doubling time scales (down to a 
few minutes) [4|. Mrk 501 has been the target of many 
multiwavelength (MWL) campaigns (e.g. Q, j6), JT), 
J8l), mainly covering the object during flaring activity. 
The data presented here were taken between March 
25th and May 16th, 2008 during an extended MWL 
campaign covering radio (Effelsberg, IRAM, Medicina, 
Metsahovi, Noto, RATAN-600, VLBA), optical (KVA), 
UV (Swift/UVOT), X-ray (RXTE/PCA, Swift/XRT and 
Swift/BAT) and 7-ray (MAGIC, Whipple, VERITAS) 
energies. The duration as well as the energy coverage 
of this particular Mrk 501 campaign are rather unique. 
Details on the participating instruments and the data 
analysis will be presented in an upcoming paper [9]. 
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Fig. 1. Combined normalized light curves for a selection of the instruments taking part in the campaign. The vertical bars denote 1 <r statistical 
uncertainties, and the horizontal bars the integration time of the observation. 



II. Light Curves 

Figure [T] shows the normalized light curve s0 far 
a selection of the instruments involved in the cam- 
paign. The different light curves cover the optical 
R band (KVA), the UV band (Swift/UVOT), the 
soft X-ray band (RXTE/PCA), the hard X-ray band 
(Swift/BAT) and the VHE band (MAGIC and VER- 
ITAS). The average fluxes for each instrument are 
given as: 4.4 mJy (KVA), 1.6 mJy (Swift/UVOT), 
2.8 • l(T 4 counts/s (Swift/BAT), 8.2 • Kr n erg/cm 2 /s 
(RXTE/PCA), 2.7 • Kr n ph/cm 2 /s (MAGIC) and 
2.2 • 10 _11 ph/cm 2 /s (VERITAS). Other instruments 
providing valuable data (like Swift/XRT or Whipple) 
have been omitted for the sake of clarity in this plot. 
Flux variations are large in X-rays and 7-rays, but rather 
small in the UV and optical. Due to the small error bars 
in the X-ray data, the most significant flux variations 
can be observed at these energies. The plot also shows 
some evidence for a correlated flux variability at X- 
rays and VHE 7-rays (see section IIVI ) indicating a 
low-emission state before MJD 54560 and a somewhat 
stronger emission afterwards. For the spectral analysis 
presented below we divided the data set into three time 
intervals taking into account the X-ray flux level (i.e. 
low/high flux before/after MJD 54560) and the data gap 
at most frequencies around MJD 54574. 

III. Variability 

We followed the description given in 1 10] to quantify 
the flux variability by means of the so-called fractional 
variability parameter F var . In order to account for the 
individual flux measurement errors (er e rr,i), the 'excess 
variance' ( ITTl . lfl2ll ') was used as an estimator of the 
intrinsic source flux variance. This is the variance after 
subtracting the contribution expected from measurement 
errors. F var was derived for each individual instrument 
taking part in the campaign, which covered an energy 

'Each individual light curve was normalized to its average. 




Fig. 2. Fractional variability parameter F var for all the instruments 
participating in the campaign. i*Var was derived using the individ- 
ual single-night flux measurements except for Swift/BAT for which 
data integrated over one week were used. Vertical bars denote 1 <r 
uncertainties, horizontal bars indicate the approximate energy range 
covered by the instrument. The arrows indicate 95% confidence level 
upper limits. 



range from radio frequencies at ~8 GHz up to very high 
energies at ~10 TeV F var is calculated as: 



IS 2 - < a 2 



> 



< F 1 > 2 



(1) 



where < F 1 > denotes the average photon flux, S 
the standard deviation of the N flux measurements and 
< <r 2 rr > the mean squared error, all determined for a 
given instrument (energy bin). The uncertainty of F val 
is estimated according to: 



AF var = 



9 < TP >2 F2 



(2) 



Fig. 12 shows the F val values derived for all instru- 
ments that participated in the MWL campaign. Some in- 
struments showed a negative excess variance (< cr 2 rr > 
larger than S 2 ), which can happen when there is little 
variability and/or the errors are slightly overestimated. 
Essentially such a result can be interpreted as no 
signature for variability in the data of that particular 
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instrument, either because a) there was no variability 
or b) the instrument was not sensitive enough to detect 
it. In these cases, upper limits of 95% confidence level 
were computed. 

The plot, on the other hand, also shows significant 
variability detected with various other instruments dur- 
ing the campaign. Essentially all instruments observing 
at optical or larger frequencies recorded variability. 
The plot also shows some evidence that the recorded 
flux variability increases with energy: in the optical R 
band (ground-based telescopes) and the 6 filters from 
Swift/UVOT the variability is around 2-4%, in X-rays it 
is about 13%, and at VHE at the 20% level, although 
affected by large error bars (due to the large uncertainties 
in the flux measurements). The radio instruments show 
no evidence for variability, with the exception of RATAN 
(22 GHz) and Metsiihovi (37 GHz) that show - 7±2%. 
In the synchrotron self-Compton (SSC) framework, the 
observed flux variability contains information on the 
dynamics of the underlying population of relativistic 
electrons (and possibly positrons). In this context, the 
general variability trend reported in Fig. [2] suggests that 
the flux variations are produced by the injection of 
energetic particles, which are characterized by shorter 
cooling time scales, causing the higher variability am- 
plitude observed at the highest energies. 

IV. MULTIFREQUENCY CROSS-CORRELATIONS 

In order to study the multifrequency cross-correlations 
between the different energy bands we used the Discrete 
Correlation Function (DCF) as described in fl3l . This 
method can also be applied in the case of unevenly 
sampled data as taken in this campaign. 
The DCF was derived for all different combinations of 
instruments / energy regions and also for artificially in- 
troduced time lags (ranging from -8 to +8 days) between 
the individual light curves. Such time lags may occur 
as a result of spatially separated emission regions of the 
individual flux components (as expected, for example, in 
external inverse Compton models), or may be caused by 
the energy dependent cooling time-scales of the emitting 
electrons. 

Based on the MWL data from this campaign, significant 
correlations have been found for the pairs RXTE/PCA 
- Swift/XRT and also (less significant) RXTE/PCA (or 
Swift/XRT) with MAGIC and VERITAS (Fig. E3 and 
[3bl . In both cases, the DCF maximum is obtained for a 
zero time lag with a value of 0.87 ±0.28 (RXTE/PCA - 
Swift/XRT ) and 0.5 ± 0.19 (RXTE/PCA - MAGIC and 
VERITAS) respectively. Due to the modest flux variabil- 
ity and / or large flux errors, no strong conclusions could 
be drawn from this analysis. 

V. SED MODELING 

The broadband SED of Mrk 501 for the three dif- 
ferent time periods defined above, together with some 
historical data from the 2005 low and high state of the 
object are shown in Fig. |4] The host galaxy contribution 
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Fig. 3. Discrete Correlation Function for time lags from -8 to +8 
days in steps of ~ 2 days. The grey band represents the expected 
fluctuation of the DCF values in the case of completely uncorrelated 
time series given the error bars from the actual observations. 

TABLE I 

THE SSC MODEL PARAMETERS USED TO DESCRIBE THE 
BROADBAND SED FOR DIFFERENT FLUX STATES OF THE 
CAMPAIGN. 





2008 


2008 


2005 


2005 




high state 


low state 


high 


low 


Tbrcak 


2.6 • 10 5 


2.2 • 10 5 


1.0 ■ 10 6 


1.0 • 10 5 


ni 


2.0 


2.0 


2.0 


2.0 


n2 


3.9 


4.2 


3.9 


3.2 


B [G] 


0.19 


0.19 


0.23 


0.31 


K [cm" 3 ] 


1.8 ■ 10 4 


1.8 ■ 10 4 


7.5 ■ 10 4 


4.3 ■ 10 4 


R [cm] 


3 ■ 10 15 


3 ■ 10 15 


1 ■ 10 15 


1 ■ 10 15 


5 


12 


12 


25 


25 



(12.0 ± 0.3 mJy [14)) has been subtracted from the 
optical (KVA) data while the 7-ray spectra have been 
corrected for EBL absorption using the Tow-IR' model 
of fl5l . The results from a one-zone SSC model fit to the 
different data sets are also shown in the figure as dashed 
lines. The model code was developed by Tavecchio et al. 
( ifTTjl . 0) and is based on the following characteristic 
parameters: a spherical emission region with radius R 
and Doppler factor 6, a magnetic field of strength B, 
an electron distribution (density K) following a broken 
power law with slopes n\ and n% and break energy 
Tbreak- The actual values of these model parameters for 
the two different emission states during the campaign 
and the historical data from 2005 are given in Tab. H] As 
can be seen from Fig. 2] the model is able to accurately 
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Fig. 4. Broadband SED for Mrk 501 as obtained during this campaign in comparison to low and high states from 2005. The results from a 
SSC model fit to the low state 2008 data (dot-dashed cuve), the high state 2008 data (heavy-dashed curve), the 2005 low state (light-dashed 
curve) and the 2005 high state (solid line curve) are shown. 



reproduce the data at X-ray energies. Given the relatively 
small differences in the SEDs of the two emission states 
of the campaign, only marginal changes of the model 
parameters were required in order to adjust the model 
to the two states. The proposed explanation for the low - 
high state transition is the injection of fresh, high-energy 
electrons which lead to a shift of the 7broak energy and 
to a hardening of the spectrum. 

The discrepancy between the model and the data at lower 
energies (radio, optical) can be caused by synchrotron 
radiation from additional, cooler electron populations 
which could be present at different locations in the 
jet. The higher (than expected) fluxes at radio/optical 
frequencies were discussed in the past (also with Mrk501 
data) in the framework of the helical-jets in blazar 
scenarios ifTTll or the blob-in-jet scenario [18|. As is 
shown in table Q] in comparison to the historical 2005 
SED, the model parameters have changed significantly. 
However, it is worth noticing that the sparse coverage of 
the 2005 data allow for a lot of degeneracy among the 
(large) number of model parameters. A robust statement 
from the comparison of the 2005 and 2008 SEDs is 
that, while the X-ray and gamma-ray fluxes did change 
substantially between these two epochs, the fluxes at 
optical frequencies remained approximately the same. 
In the framework of two populations of electrons, this 
result suggests that the population of cool electrons does 
not vary with time while the population of electrons 
responsible for the X-ray (Synchrotron) and gamma- 
ray (Inverse Compton) emission is very dynamic. A 
more detailed modeling of the experimental data will 
be performed in a forthcoming publication |9|. 
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